Matrix Extension with Symmetry and Construction of 
Biorthogonal Multi wavelets 



Xiaosheng Zhuang*'^ 

"Department of Mathematical and Statistical Sciences, University of Alberta, Edmonton, Alberta, Canada 

T6G 2G1. 



Abstract 

Let (P, P) be a pair of r x i matrices of Laurent polynomials with symmetry such that 
P(z)P*(z) = Ir for all z e C\{0) and both P and P have the same symmetry pattern 
that is compatible. The biorthogonal matrix extension problem with symmetry is to 
find a pair of s x s square matrices (P^,, P^.) of Laurent polynomials with symmetry 
such that [/r, 0]P;. - P and [7^,0] Pg = P (that is, the submatrix of the first r rows 
of Pf,, Pe is the given matrix P, P, respectively), P^ and P^, are biorthogonal satisfying 
Pj,(z)P*(z) - Is for all z € C\{0}, and P^,, Pe have the same compatible symmetry. 
In this paper, we satisfactorily solve this matrix extension problem with symmetry 
by constructing the desired pair of extension matrices (P^, P^) from the given pair of 
matrices (P, P). Matrix extension plays an important role in many areas such as wavelet 
analysis, electronic engineering, system sciences, and so on. As an application of our 
general results on matrix extension with symmetry, we obtain a satisfactory algorithm 
for constructing symmetric biorthogonal multiwavelets by deriving high-pass filters 
with symmetry from any given pair of biorthgonal low-pass filters with symmetry. 
Several examples of symmetric biorthogonal multiwavelets are provided to illustrate 
the results in this paper 
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1. Introduction and Main Result 

The matrix extension problem plays a fundamental role in many areas such as elec- 
tronic engineering, system sciences, mathematics, and etc. To mention only a few 
references here on this topic, see 



example, matrix extension is an indispensable tool in the design of filter banks in elec- 
tronic engineering ( LI61, 2^, 2l]]) and in the construction of multiwavelets in wavelet 
analysis (H H S larTTli l'l li Il5l [itI IH]). In order to state the biorthogonal matrix 
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extension problem and our main result on this topic, let us introduce some notation and 
definitions first. 

Let p(z) = 2ifcez Pkz'', z e C\{0) be a Laurent polynomial with complex coefficients 
Pk & C for all k e Z. We say that p has symmetry if its coefficient sequence {pk]kez has 
symmetry; more precisely, there exist e e {-1, 1} and c e Z such that 

Pc-k^spk, yke-L. (LI) 

If e = 1, then p is symmetric about the point c/2; if e = -1, then p is antisymmetric 
about the point c/2. Symmetry of a Laurent polynomial can be conveniently expressed 
using a symmetry operator S defined by 

^pa):=-7j^, zeC\{0). (L2) 
P(l/z) 

When p is not identically zero, it is evident that ( ILlb holds if and only if Spiz) = 
sz'^- For the zero polynomial, it is very natural that SO can be assigned any symmetry 
pattern; that is, for every occurrence of SO appearing in an identity in this paper, SO 
is understood to take an appropriate choice of sz'^ for some s e {-1,1} and c e Z so 
that the identity holds. If P is an r x s matrix of Laurent polynomials with symmetry, 
then we can apply the operator S to each entry of P, that is, SP is an r x s matrix 
such that [SP]j^k := >S([P]j;t), where [P]ik is the (j,k)-entry of the matrix P. Also 
[P]j.k:e := [[P]j,k, [P]jMu . . . , [P]j,e] isalx({-k+l) vector 

For two matrices P and Q of Laurent polynomials with symmetry, even though all 
the entries in P and Q have symmetry, their sum P + Q, difference P - Q, or product 
PQ, if well defined, generally may not have symmetry any more. This is one of the 
difficulties for matrix extension with symmetry. In order for P + Q or PQ to possess 
some symmetry, the symmetry patterns of P and Q should be compatible. For example, 
if SP - SQ, that is, both P and Q have the same symmetry pattern, then indeed 
P + Q has symmetry and S{P + Q) = SP = SQ. In the following, we discuss the 
compatibility of symmetry patterns of matrices of Laurent polynomials. For an r x i 
matrix P(z) - Yjk^z Pkz'', we denote 

P'-(z):=J^p;z-'' with P;:^pf, keZ, (1.3) 

keZ 

T 

where Pk denotes the transpose of the complex conjugate of the constant matrix in 
C. We say that the symmetry ofP is compatible or P has compatible symmetry, if 

SP(z) = (Seinz)S02(z), (L4) 

for some 1 x r and 1 x .? row vectors Oi and 02 of Laurent polynomials with symmetry. 
For an rx s matrix P and an i x f matrix Q of Laurent polynomials, we say that (P, Q) 
has mutually compatible symmetry if 

SP(z) ^ (Seinz)Se(z) and SQ(z) ^ iSenz)Se2(z) (1.5) 

for some Ixr, Ixi, Ixf row vectors 6i , 0, 02 of Laurent polynomials with symmetry. 
If (P, Q) has mutually compatible symmetry as in (IL5b . then their product PQ has 
compatible symmetry and in fact >S(PQ) = (S0iyS62. 
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For a matrix of Laurent polynomials, another important property is the support of 
its coefficient sequence. For P = TjkezPkz'' such that Pk - for all k e 1\{m,n\ with 
P,n + and P„ + 0, we define its coefficient support to be csupp(P) :— [m,n] and 
the length of its coefficient support to be |csupp(P)| :- n - m. In particular, we define 
csupp(O) :- 0, the empty set, and |csupp(0)| := -oo. Also, we use coeff(P,fe) := Pk to 
denote the coefficient matrix (vector) Pk of in P. Throughout this paper, always 
denotes a general zero matrix whose size can be determined in the context. 1„ denotes 
the 1 X « row vector [1, . . . , 1], 

The Laurent polynomials that we shall consider have their coefficients in a subfield 
F of the complex field C. Several particular examples of such subfields F are F = Q 
(the field of rational numbers), F = R (the field of real numbers), and F = C (the field 
of complex numbers). 

Throughout the paper, r and s denote two positive integers such that 1 ^ r ^ i. 
Now we generalize the matrix extension problem we consider in Iil4ll to the biorthog- 
onal case as follows: Let (P, P) be a pair of r y. s matrices of Laurent polynomials 
with coefficients in F such that P(z)P'(z) - Ir for all z e C\{0}, the symmetry of 
each P and P is compatible, and SP - SP. Find a pair of s x s square matrices 
(Pf,, Pe) of Laurent polynomials with coefficients in F and with symmetry such that 
[/r,0]P(, - P, [/r,0]Pf. = P, (that is, the submatrix of the first r rows of P,., P^ is 
the given matrix P, P, respectively), the symmetry of Pe and Pe is compatible, and 
Pj,(z)P*(z) - Is for all z e C\{0}. The coefficient support of Pe, Pe can be controlled by 
that of P, P in some way. 

The above extension problem plays a critical role in wavelet analysis. The key of 
wavelet constructions is the so-called multiresolution analysis (MRA), which contains 
mainly two parts. One is on the construction of refinable function vectors that satisfies 
certain desired conditions. For example, (bi)orthogonality, symmetry, regularity, and 
so on. Another part is on the derivation of wavelet generators from refinable function 
vectors obtained in first part, which should be able to inherit certain properties similar 
to their refinable function vectors. From the point of view of filter banks, the first part 
corresponds to the design of filters or filter banks with certain desired properties, while 
the second part can be and is formulated as a matrix extension problem given above. 
For the construction of biorthogonal refinable function vectors (a pair of biorthogo- 
nal low -pass filters), the CBC {coset by coset) algorithm proposed in |10] (also see 
Section 3 for more details) provides a systematic way of constructing a desirable dual 
mask from a given primal mask that satisfies certain conditions. More precisely, given 
a mask (low-pass filter) satisfying the condition that a dual mask exists, following the 
CBC algorithm, one can construct a dual mask with any preassigned orders of sum 
rules, which is closely related to the regularity of the refinable function vectors. Fur- 
thermore, if the primal mask has symmetry, then the CBC algorithm also guarantees 
that the dual mask has symmetry. Thus, the first part of MRA corresponding to the 
construction of biorthogonal multiwavelets is more or less solved. However, how to 
derive the wavelet generators (high-pass filters) with symmetry remains open even for 
the scalar case {r = 1) and this is one of the motivations of this paper We shall see 
that using our extension algorithm, the wavelet generators do have symmetry once the 
given refinable function vectors possess certain symmetry patterns. 
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Due to the flexibility of biorthogonality PP* - Ir, the above extension problem 
becomes far more complicated than that the matrix extension problem we considered 
in II14II . The difliculty here is not the symmetry patterns of the extension matrices, but 
the support control of the extension matrices. Without considering any issue on sup- 
port control, almost all results of Theorems 1 and 2 in fTJ] can be transferred to the 
biorthogonal case without much difficulty. In lll4ll . we showed that the length of the 
coefficient support of the extension matrix can never exceed the length of the coeffi- 
cient support of the given matrix. Yet, for the extension matrices in the biorthogonal 
extension case, we can no longer expect such nice result, that is, in this case, the length 
of the coefficient supports of the extension matrices might not be controlled by one of 
the given matrices. Let us present an example here to show why we might not have 
such a result. 



Example 1. Consider two 1x3 vectors of Laurent polynomials p(z) - [1,0, a{z)] and 
p(z)^ [l,a'(z),0] with |csupp(a(z))| > 0, |csupp(fl(z))| > 0. We have_pp* = 1. Let P^ 
and Pf be their extension matrices such that P(,P* = 73. Then P^., Pg must be of the 
form: 



P. = 



1 

-b2(z)a*(z) 





bdz) 
biiz) 



a(z) 
ciiz) 
ciiz) 



1 a{z) 

-ci(z)fl*(z) biiz) c'i(z) 

-C2(z)fl*(z) biiz) C2(z) 



It is easy to show that det(P£.) - b\{z)c2{z) - b2(z)ci(z). Since P^ is invertible with 
P^' - P*, we know that det(P(,) must be a monomial. Without loss of generality, we 
can assume b\{z)c2{z) - b2{z)ci{z) - I. Using the cofactors of P^, it is easy to show 
that Pe - (Pj')' must be of the form: 



P. = 



1 a(z) 

bl{z)a*(z) ciiz) + a(z)a''(z)b]{z) -b^iz) 
-bl{z)a*(z) -cl(z) -a(z)a*(z)bl(z) b^iz) 



On the one hand, if |csupp(Z7i(z))| > or |csupp(Z72(z))| > 0, then we see that one 
of the extension matrices will have support length exceeding the maximal length of the 
given columns. One the other hand, if both |csupp(Z7i(z))| - and |csupp(/72(z))l = 
(in this case, both bi(z) and b2{z) are monomials), then the lengths of the coefficient 
support of ci(z) and C2(z) in Pg must be comparable with 'a*(z)a{z) so that the support 
length of Pe can be controlled by that of p or p, which in turn will result in longer 
support length of Pg. 



The above example shows that it is difficult to control the support length of the 
coefficient support of the extension matrices independently by only one given vector 
in the biorthogonal setting. Nevertheless, we have the following result, which indicate 
the lengths of the coefficient support of the extension matrices can be controlled by the 
given pair in certain sense. 

Theorem 1. Let ¥ be a subfield ofC Let (P, P) be a pair of r x s matrices of Laurent 
polynomials with coefficients in F such that the symmetry of each P, P is compatible: 
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SP — SP — (S6i yS62 for some 1 X r, 1 X vectors 6i , 62 of Laurent polynomials with 
symmetry. Moreover, P(z)P*(z) = lyforall z e C\{0). Then there exists a pair of s X s 
square matrices (Pg, Pj) of Laurent polynomials with coefficients in F such that 

(i) [7^,0] Pf — P, [Ir, 0]Pe — P, that is, the submatrices of the first r rows ofPg, Pg 
are P, P, respectively; 

(ii) Pe and P^ are biorthogonal: Pe(z)Pl(z) — Is for all z e C\{0}," 

(iii) The symmetry of each Pe,Pe is compatible: SPe — SPe — (SO)* SO2 for some 
\ s vector 6 of Laurent polynomials with symmetry. 

(iv) Pf,, Pe can be represented as: 

Pe(z)^Pj{z)---Px{z), P,(Z)= Py(z)---Pi(z), (1.6) 

where P j, P j, 1 ^ 7 ^ 7 are sxs matrices of Laurent polynomials with symmetry 
that satisfy Pj(z)P*(z) = /j. Moreover, each pair of(Pj+i, Pj) and (Pj+i, Pj) has 
mutually compatible symmetry for all j — I, . . . , J — I. 

(v) If r — 1, then the coefficient supports of P^, Pe <^re controlled by that of P, P in 
the following sense: 

max {|csupp([P,];.j)|,|csupp([P,]j.t)|) max |csupp([P]f)| + max |csupp([P],)|. 

(1.7) 

For r = 1, Goh et al. in |8] considered this matrix extension problem without 
symmetry. They provided a step-by-step algorithm for deriving the extension matrices, 
yet they did not concern about the support control of the extension matrices nor the 
symmetry patterns of the extension matrices. For r > 1, there are only a few results 
in the literature Jll |4|] and most of them concern only about some very special cases. 
The difficulty still comes from the flexibility of the biorthogonality relation between 
the given two matrices. In this paper, we shall mainly consider this matrix extension 
problem with symmetry for the biorthogonal case and shall provide an extension algo- 
rithm from which the extension matrices can have both symmetry and support control 
as stated in Theorem[T] 

Here is the structure of this paper In Section 2, we shall introduce some auxiliary 
results, prove Theorem[T] and also provide a step-by-step algorithm for the construction 
of the extension matrices. In Section 3, we shall discuss the applications of our main 
result to the construction of symmetric biorthogonal multiwavelets in wavelet analysis. 
Examples will be provided to illustrate our algorithms. Conclusions and remarks shall 
be given in the last section. 

2. Proof of Theorein[T]and an Algorithm 

In this section, we shall prove our main result Theorem [T] and based on the the 
proof, we shall provide a step-by-step extension algorithm for deriving the desired pair 
of extesion matrices. 

First, let us introduce some auxiliary results. The following lemma shows that 
for a pair of constant vector (f , f) in F, we can find a pair of biorthogonal matrices 
(^(f f)' ^(f f)) ^^'-^ '^hat up to a constant multiplication, they normaUze £, £ to two unit 
vectors, respectively. 
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Lemma 1. Let (f , f) be a pair of nonzero 1 x n vectors in F. Then the following 
statements hold. 

(1) If ff* + 0, then there exists a pair of n X n matrices (u^f^y ^(f f)) ^ ^uch 
that = [(iy,F], [/(f = [(f)%F], ant/ t/(f,f)C/*f^ = In, where F,F are 
n X (n — 1) constant matrices in ¥ and c,~c are two nonzero numbers in F such 
that f f* = cc. In this case, £U^^j^ — cCi and f f^(f f) = cCi. 

(2) If ff* = 0, then there exists a pair of n X n matrices {^U^^jy ^(f f)) ^ such 

thatU^,j^j U^,~,^ = a««f t/(,_^,f7*^_ = 4, 

where F, F are n X (n — 2) constant matrices in F and ci, C2, c^i, Q are nonzero 
numbers in F iMc/; f/zaf ||f|p - c{ci,\\f\\^ — C2C2. /n this case, £t^(f f) = ciCi and 

^^(f,f) = '^262- 

Proof. If f £* 0, there exists {£2, . . . , f,,} being a basis of the orthogonal compliment 
of the linear span of {£) inF". LetF:= [£;,...,£*] andU^^j^ := [(f)*, F]. Then t/^^ 

is invertible. Let [7,^7, := . It is easy to show that t/(f fj and t/(f fj are the 

desired jnatrices. 

If ££* = 0, let {£3, ... , £„) be a basis of the orthogonal compliment of the linear 
span of {£, £} in F«. Let U^^j^ = [(|)% with F [£;, . . . , £^1. Then U^,j^ 

andt/,^T, := are the desired matrices. ■ 

Thanks to Lemma [T] we can reduce the support lengths of a pair (p, p) of Laurent 
polynomials with symmetry by constructing a pair of biorthogonal matrices (B, B) of 
Laurent polynomials with symmetry as stated in the following lemma. 

Lemma 2. Let (p, p) be a pair oflxs vectors of Laurent polynomials with symmetry 
such that pp* - 1 and Sp - Sp - ez'^[ls,, -1.52,z"'lj3, -z^'ljj =: SO for some 
nonnegative integers ii, . . . , 54 satisfying si + ■ ■ ■ + S4 — s and e € {1, —l],c € {0, 1). 
Suppose |csupp(p)| > 0. Then there exist a pair of s X s matrices (B, B) of Laurent 
polynomials with symmetry such that 

(1) B,B are biorthogonal: B(z)B*(z) = I„; 

(2) .SB = .SB = (SOySOi withSOi = szf'[l,,'^,-l,,'^,z~^ls'^,-z~^ls'J for some non- 
negative integers s'^, . . ., s'^ such that s'^ + ■ ■ ■ + s'^ — s; 

(3) the length of the coefficient support of p is reduced by that o/B. B does not in- 
crease the length of the coefficient support of p. That is, |csupp(pB)| ^ |csupp(p)|- 
|csupp(B)| and |csupp(pB)| ^ |csupp(p)|. 

Proof. We shall only prove the case that SO = [1^,, -lj,,z"^lj3, -z^^ls^]- The proofs 
for other cases are similar By their symmetry patterns, p and p must take the form as 
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t-1 



follows with > and coeff(p, -€) + 0: 

p = [fi,-£2,gi,-g2k"^ + [£3,-£4,g3,-g4k"^^' + ^ coeflF(p,^)/ 
+ [£3,£4,gi,g2]/"' +[fi,f2,0,0k^; 

_____ - _ '?-2 

p = [£i, -f2, gi, -g2]z"'' + [f3, -f4, g3, -g4]z"''^' + ^ coeff(p,^)z* 
+ [f3,f4,gi,g2]z'"' + [fi,f2,0,0k'; 



(2.1) 



Then, either ||£i|| + ||£2|| or ||gi|| + HgjH 9^ O^Considering ||£i|| +J£2|| 9^0, due to 
pp* = 1 and |csupp(p)| > 0, we have £i£* - £2£2 - 0. Let c :- £i£* = £2£2. Then 
there are at most three cases: (a) c 0; (b) c = but both £i, £2 are nonzero vectors; 
(c) c = and one of £1, £2 is 0. __ __ 

Case (a): In this case, we have £:£* and £2 £2 + 0. By Lemma[T] we can con- 
struct two pairs of biorthogonal matrices (f^(f| f^j, ^^(f, f,)) ^i^^ (f^(f, f^j, ^(fzfj)) ^^'-^ 
respect to the pairs (£1, £1) and (£2, £2) such that 



^(f2,f2) ~ 









Si, 






±2 






Si, 


,F2 





(fi.fi) 



where c\ , ci are constants in F such that c - c\C\. Define Bo(z), Bo(z) as follows: 
Bo(z) 



Bo(z) = 



l+J ' / f] \* 
2 "-eL-^ 


^1 


2 "-^i 








2 





2 ^7, 


Fi 





[ 













2 V, 


Fx 


2 ^c, 








2 Vi 





2 c, 


Fi 





[ 











/.S, + .!4 _ 



(2.2) 



Direct computation shows that Bo(z)B()(z)* - I, due to the special structures of the 
pairs j^^) and (t/(fj f,), t^(f, f^)) constructed by Lemma[T] The symmetry 

patterns of pBo and pBo satisfies 



.S(pBo) = .S(pBo) = [z-',l„-i,-z- 



l,,_l,Z-4,3,-Z-'l,J. 



Moreover, Bo(z), Bq(z) reduce the lengths of the coeflicient support of p, p by 1, re- 
spectively. 

In fact, due to the above symmetry pattern and the structures of Bq, Bq, we only 
need to show that coeff'([pBo];, ^) - coeff'([pBo] ^) = for 7 = 1, ii -1- 1. Note 
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that coeff([pBo]/,^) = coeff(p,^)coeff([Bo]. i,0) = i(fif! - fzf;) = 0. Similar 

computations apply for other terms. Thus, |csupp(pBo)| < csupp(p) and |csupp(pBo)| < 
|csupp(p)|. Let £■ be a permutation matrix such that 



^(pBo)£ = 5(pBo)£ = 



.S3+1, 



S4+1 



] =: se. 



Define B(z) - BqCz)/? and B(z)^ = Bo(z)£'. Then B(z) and B(z) are the desired matrices. 

Case (b): In this case, £if* = £2£5 =0 and both £1, £2 are nonzero vectors. We 
have £i£j + and £2 £2 + 0. Again, by Lemma [T] we can construct two pairs of 
biorthogonal matrices (t/(f,,f,), t/(f,,fi)) and ({/(f^.fj), t/(f2,f2)) with respect to the pairs 
(£1, £1) and (£2, £2) such that 



f^(fi,fi) = 

f^(f2,f2) = 



Cl 



' t^(fl,fl) - 

' f^(f2,f2) = 



Co 
Co/ 



, £it/(fi,fi) = coci, 

, £2f/(f,,f2) = CoCi, 



where co, ci, C2 are constants in F such that £i£j = coci and £2£2 = coC2. Let Bo, Bo(z) 
be defined as follows: 





2 '-^i > 


Fx 


2 








Bo(z) = 


2 "-o 







Fi 






















2 )ca> 


i^i 


2 '-Co'' 








Bo(z) = 


2 l^co^ 





2 "-en-* 


Fi 






















(2.3) 



We can show that Bo(z) reduces the length of the coefficient support of p by 1, while 
Bo(z) does not increase the support length of p. Moreover, similar to case (a), we can 
find a permutation matrix E such that 



S(pBo)E = S(pBo)E = [1 



1,-1, 



S3+1, 



S4+1 



] =: S6i 



Define B(z) - Bo(z)£' and B(z) = Bo{z)E. Then B(z) and B(z) are the desired matrices. 

Case (c): In this case, £i£* = £2£2 = and one of £1 and £2 is nonzero. Without 
loss of generality, we assume that £1 and £2 = 0. Construct a pair of matrices 
(^(fi fi)' ^(fi fi)) Lemma[T]such that fit^(f| f^j - ciCi and £it^(f| f^j - C2t2 (when 
£1 = 0, the pair of matrices is given by (t/(f,,f|), t/(f,,f,)))- Extend this pair to a pair of 
sxs matrices {U, U) by U := diag (f/(f,_f,), and U diag (f/(f,_f,), Then 
pf/ and pt/ must be of the form: 



q := pt/ = [ci,o,...,0,-£2,gi,-g2k + [£3,-£4,g3,-g4jz 

C-2 

+ coeff(q,^)z^ + [£3,£4,gi,g2]/"' +[ci,0,...,0,£2,0,0]z'; 



k=-l+2 
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q pf/ = [0, C2, . . . , 0, -£2, gi, -gilz"' + [£3, -U, 93, -94]^"'^' 

7-2 

+ ^ coeflF(q, + [£3, £4, 91, g2]z'"' + [0, cj, . . . , 0, , £2, 0, 0]z'; 

k=-l+i 

If = 0, we choose k such that k - argminf^i{|csupp([q]i)| - |csupp([q]f)|), i.e., k 
is an integer such that the length of coefficient support of |csupp([q]i)| - |csupp([q]j^)| 
is minimal among those of all |csupp([q]i)| - |csupp([q]f)|, { - 2, . . . , s; otherwise, due 
to qq* - 0, there must exist k such that 

|csupp([q]i)| - |csupp([q]i)K max |csupp([q];)| - |csupp([q]i)|, 

(k might not be unique, we can choose one of such k so that |csupp([q] i)| - |csupp([q]jt)| 
is minimal among all |csupp([q]i)| - |csupp([q]f)|, { = 2, . . . , s). 

For such k (in the case of either [q] 1 =0 or [q] 1 0), define two matrices B(z), B(z) 
as follows: 



1 


• 


• 






1 


■■ 


• a(zy 







1 ■ 


■ 




,B(z) = 





1 ■■ 







-a(z) 


■ 


■ 1 









■■ 


1 












Is-k 



where a(z) in B(z), B(z) is a Laurent polynomial with symmetry such that Sa(z) — 
S{{ai\i)ISi{ai\k), |cmpp([q]i -fl(z)[q]t)| < |csupp([q]^)|, and |csupp([q]t -a(z)*[q]i)| < 
maxi^f^j |csupp([q]f)|. Such a{z} can be easily obtained by long division. 

It is straightforward to show that B(z)B*(z) = B(z) reduces the length of the 
coefficient support of q by that of a{z) due to |csupp([q]i - fl(z)[q]t)| < |csupp([q]i:)|. 
And by our choice of k, B(z) does not increase the length of the coefficient support of 
q. Moreover, the symmetry patterns of both q and q are preserved. 

In summary, for all cases (a), (b), and (c), we can always find a pair of biorthogonal 
matrices (B, B) of Laurent polynomials such that B reduces the length of the coefficient 
support of p while B does not increase the length of the coefficient support of p. 

For ||£i|| + ||£2|| - 0, we must have ||gi|| + ||g2|| 0. The discussion for this case is 
similar to above. We can find two matrices B(z), B(z) such that all items in the lemma 
hold. In the case that gigf = g2g2 = c\c\ + 0, the pair (Bo(z), Bo(z)) similar to ( 12.2b is 
of the form: 



















Bo(z) = 





1+z / 91 \* 
2 ^■'SyJ 


Gi 


i-z/ai -1* 
2 > 










2 "-^i > 





i+z/'g2\* 
2 "-^i > 


G2 




^.!l+.!2 











1 


Bo(z) = 





1+z/Si \* 
2 Vi 


Gi 


l-z/ 3i \* 










1-2/92 \* 

2 





i+z/'g2\* 


G2 



(2.4) 



The pairs for other cases can be obtained similarly. We are done. ■ 
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Let be a 1 X « row vector of Laurent polynomials with symmetry such that SO = 
[s-iz'^' , . . . , e„z'^"] for some si,. . .,s„ e (-1,1) and ci, . . . , c„ e Z. Then, the symmetry 
of any entry in the vector 0diag{z^^'''^^\ . . . ,2"^"^"^'^) belongs to {+1, ±z"')- Thus, there 
is a permutation matrix Eg to regroup these four types of symmetries together so that 



S(e\Js9) = [i„,,-i„„z-'i„„-z-'i„j. 



,-rc„/2i 



)Eg and rii 



(2.5) 

n4 are nonnegative integers 



where Use ■- diag(z"'^'''^^l, . . . ,z 
uniquely determined by SO. 

For an rx matrix P of Laurent polynomials with compatible symmetry as in ( |1.4| i, 
it is easy to see that Q := U^g^ PU^e^ has the symmetry pattern as follows. 

SQ = [l„,-l,„zl,3,-zl,J^[l„,-l,„z-'l,3,-z-il,J. (2.6) 

Note that U^e, and U^e, do not increase the length of the coefficient support of P. 
Now, we can prove Theorem[T]using Lemma|2] 

Proof (Proof of Theorem [T]). First, we normalize the symmetry patterns of P and P 
to the standard form as in ( 12.6b . Let Q := U^g PU^e, and Q := U^g PU^e, (given 0, 

Use is obtained by ( 12.51 )). Then the symmetry of each row of Q or Q is of the form 
ez''[ls,,-lsj,z"'ls3,-z"'ljj for_someeG {-l,l}andc€ {0^1]. 

Let p :- [Q]i,: and p := [Q]i,: be the first row of Q, Q, respectively. Applying 
Lemma|2]recursively, we can find pairs of biorthogonal matrices of Laurent polynomi- 
als (Bi,Bi), {Bk,Bk) such that pBi ■■■B/f = [1,0,..., 0] andpBi ■■ - B/f = [l,q(z)] 
for some 1 x (s - 1) vector of Laurent polynomials with symmetry. Note that by 
Lemma|2l all pairs (Bj, Bj+i) and (Bj, Bj+i) for j = I, . . . , K - 1 have mutually com- 
patible symmetry. Now construct Bk+i(z), Bx+iiz) as follows: 



B/r+i(z) = 



1 

q*(z) 







Bk^i(z) = 



1 -q(z) 




Bk+\ and B^r+i are biorthogonal. Let A : 

Then pA = pA = d. 

Note that QA and QA are of the forms: 



QA = 





Qi(z) 



Bi ■ ■ ■ BkBk+1 and A := Bi ■ ■ ■ BkB 



QA = 





Qi(z) 



for some (r- 1) x s matrices Qi, Qi of Laurent polynomials with symmetry. Moreover, 
due to Lemma |2] the symmetry patterns of Qi and Qi are compatible and satisfies 
SQi = SQi. The rest of the proof is completed by employing the standard procedure 
of induction. ■ 

According to the proof of Theorem [l] we have an extension algorithm for Theo- 
rem[T| See Algorithm 1 . 
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Algorithm 1 Biorthogonal Matrix Extension with Symmetrty 

(a) Input: P, P as in Theorem[T]with SP = SP = (SeiyS62 for two 1 x r, 1 x i row 
vectors 9i , 02 of Laurant polynomials witfi symmetry. _ 

(b) Initialization: Let Q := PU^e, and Q := Uj^g^PU^e^- Then both Q and Q have 



the the same symmetry pattern as follows: 

SQ=SQ^ [lr„-lr„zlr,,-zlrf[K,-h„Z-'l,,,-Z-'h,l 



(2.7) 



where all nonnegative integers ri, . . . , r4, si, . . . , S4 are uniquely determined by SP. 

Note that this step does not increase the lengths of the coefficient support of both 

P and P. 
(c) Support Reduction: _ 
1: Let Uo := U^g and A - A :- 1^. 
2: for ^ = 1 to r do _ 
3: Let p := [Q]k,t.s and p := [Q]i,fc,^ 
4: while |csupp(p)| > and |csupp(p)| > do 

5: Construct a pair of biorthogonal matrices {B{z), B(z)j with respect to the pair 
(p, p) by Lemma|2]such that 

|csupp(pB)| + |csupp(pB)| < |csupp(p)| + |csupp(p)|. 

6: Replace p, p by pB, pB, respectively. 
7: Set A := Adiag(4_i, B) and A := Adiag(4_i, B). 
end while 

The pair (p, p) is of the form: ([1,0,..., 0], [l^q(z)]) for some lx(s -k) vector 
of Laurent polynomials q(z). Construct B(z), B(z) as follows: 



B(z) = 



1 

q*(z) 





Is-k 



,B(z) = 



1 -q(z) 
L-k 



10: Set A := Adiag(4_i, B)^nd A := Adiag(4_i, B). 
11: SetQ := QAandQ := QA. 

12: end for _ _ 

(d) Finalization: Let Ui := diagiUsg^J,^,.)- Set P, := UiA*Uo and P, := UiA*Uo. 

(e) Output: A pair of desired matrices (P^,, P^) satisfying all the properties in Theo- 
rem [T] 



3. Application to Biorthogonal Multiwavelets with Symmetry 

In this section, we shall discuss the connection between matrix extension and 
biorthogonal multiwavelets. We shall also discuss the application of our results ob- 
tained in previous section to the construction of biorthogonal multiwavelets with sym- 
metry. Several examples are provided to demonstrate our results. 

We say that d is a dilation factor if d is an integer with |d| > 1. Throughout this 
section, d denotes a dilation factor For simplicity of presentation, we further assume 
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that d is positive, while multiwavelets and filter banks with a negative dilation factor 
can be handled similarly by a slight modification of the statements in this paper. 

Let F be a subfield of C. A low-pass filter aq : Z i-> F''^'' with multiplicity r is a 
finitely supported sequence of r x r matrices on Z. The symbol of the filter oq is defined 
to be ao(z) := Yjkez cia{k)z'^, which is a matrix of Laurent polynomials with coefficients 
in F. Let d be a dilation factor and d\,d2 be two fixed number in F such that d, - d\d2 
(for instance d\ — l,d2 - 2 for d = 2 if F = Q). Let (fl(),'flo) be a pair of low-pass 
filters with multiplicity r. We say that (oq, oo) is a pair of biorthogonal di-band filters if 



d-l 

^ ao;y(z)ao.y(z) = 

r=o 



z e C\{0}, 



(3.1) 



where ao;y and ao;y are d-band subsymbols (polyphases, cosets) of ao and ao defined to 
be 

ao;y(z) := di Zkezooik + dk)z'', ^ ^ ,^ 
ao-y{z) -.^ diY,kezao{k + dk)z^, 

Quite often, a low-pass filter aq is obtained beforehand. To construction a pair of 
biorthogonal d-band filters (aq, «()), i.e., (13.1b holds, one can use the CBC {coset-by- 
coset) algorithm proposed in [ 10] to derive tio from aq. There are two key ingredients 
in the CBC algorithm. One is that the CBC algorithm reduces the nonlinear system in 
the definition of sum rules for Hq to a system of linear equations. Another is that the 
CBC algoirithm reduces the big system of linear equation of biorthogonality relation 
for the pair (aq, a^o) to a small systems of hnear equations in (13. 11 1. Moreover, the CBC 
algorithm guarantees that for any given positive integers^, there always exists a finitely 
supported filter a^o that satisfies the sum rules of order ^. For more details on the CBC 
algorithm, one may refer to ifioi [12I1 . In our example presented in this section, the 
pairs of biorthogonal d-band low-pass filters are obtained via this way using the CBC 
algorithm (see examples in [T?]). 

For / e Li(M), the Fourier transform used is defined to be /(^) := j^f{x)e^"^^dx 
and can be naturally extended to L2(M) functions. For a pair of biorthogonal d-band 
filter (ao,a()), we assume that there exist a pair of biorthogonal d-refinable function 
vectors {<p, 4>) associated with the pair of biorthogonal d-band filters (ao, ao). That is. 



^*(d^) = ao(e-'^)0(a, <^(d^) = ade-'^)m ^ e 



and 



{,p(--k),,p):-- 



(x - k)^(x) dx = 6(k)I,; 



keZ, 



(3.3) 



(3.4) 



where 6 denotes the Dirac sequence such that (5(0) = 1 and 6{k) = for all k i^Q. 

To construct biorthogonal multiwavelets in L2(]R), we need to design high-pass 
filters Ai, . . . , Ad-i : Z F''^'' and Hi, . . . , Ad-i '■ Z — > F''^'' such that the polyphase 
matrices with respect to the filter banks {ao, ai , . . . , aj-i ) and fao, ai , . . . , aj-i ) 



P(z) = 



ao;o(z) 
ai;o(z) 

ad-i;o(z) 



ao;d-i(z) 
ai;d-i(z) 



ad-l;d-l(z) 



P(z) = 



ao;o(z) 
ai;o(z) 



ad-i;o(z) 



ao;d-i(z) 
ai;d-i(z) 



ad-i;d-i(z) 



(3.5) 
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are biorthogonal, that is, P(z)P*(z) - 7^^, where a„,-y, a,„;y are subsymbols of a,„, 
defined similar to (13.2b for m,j - 0, . . . ,d - 1, respectively. The pair of filter banks 
({ao, . . . , ad-i ), [ao, . . . , a^-i )) satisfying PP* - 7^^ is called a pair of biorthogonal 
filter banks with the perfect reconstruction property. 

Symmetry of the filters in a filter bank is a very much desirable property in many 
applications. We say that the low-pass filter ao (or ao) has symmetry if 

ao(z) = diag(eiz^''' , . . . , e,z^^Oao(l /z)diag(eiz-^' , . . . , (3.6) 

for some e\,. . .,Er e {-1,1} and ci , . . . , c,. e M such that dcf - cy € Z for all (, j - 
1 , . . . , r. If ao has symmetry as in (13. 6t and if 1 is a simple eigenvalue of ao( 1 ), then it is 
well known that the d-refinable function vector (p in (I3.3l l associated with the low-pass 
filter ao has the following symmetry: 

(^l(ci - ■) = £101, 02(C2 - •) = £2(^2, ^r(Cr - ■) ^ Sr4>r- (3.7) 

Under the symmetry condition in (|3.6l l, to apply Theorem[Tl we first show that there 
exists a suitable invertible matrix U, i.e., det(U) is a monomial, of Laurent polynomials 
in F acting on := [ao;o, . . . , ao;d-i] so that PaoU has compatible symmetry. Note 
that Pa„ itself may not have compatible symmetry. 

Leiiuna 3. Let Pa^ := [ao;o, • . ■ , ao;d-i], where ao;o, ■ • . , ao;ci-i are d-band subsymbols 
of a low-pass filter ao satisfying (I3.6l l. Then there exists a dr xdr invertible matrix U 
of Laurent polynomials with symmetry such that PaoU has compatible symmetry. 

Proof. From ( I3.6l l, we deduce that 

[ao;r(z)]c; = e^s/'' [aa-Qliz-')]t.j, r = 0, . . . , d - 1 ; ; = 1 , . . . , r, (3.8) 
where . eY :- [Q, . . . ,d - I] and R^^ ., Q^^ . are uniquely determined by 

dee -Cj-y^ dR]. + Q\. with R\. e Z, Q]. e T. (3.9) 

Since dee - cj e Z for all {,j - 1 , . . . , r, we have ce - cj e Z for all ^, y = 1 , . . . , r 
and therefore, Q^^ j is independent of £. Consequently, by (13.81 1. for every I ^ j ^ r, 
the jth column of the matrix ao;y is a flipped version of the jth column of the matrix 
ao;2i: . Let Kj^y e Z be an integer such that |csupp([ao;y]:j +z'^''''[ao;Q>: .]:j)| is as small as 
possible. Define P :- [bo;o, . . . , bo;d-i] as follows: 



[bo;y]:J : = 



5([a0;r]:,> + Z^-[ao;eJ,]:j), r < Gf,/ (3.10) 



i([aO;y]:,> - z''^' [ao;er ,],y), y > Q] ., 



where [ao;y]:,; denotes the jth column of ao;y. Let U denote the unique transform matrix 
corresponding to ( 13.101 ) such that P := [bo;o, ■ . . , bo;d-i] - [ao;o, . . . , ao;d-i]U. It is 
evident that U is paraunitary and P - PaoU. We now show that P has compatible 
symmetry. Indeed, by ( I3.8l l and ( 13.10b . 

[SbQ;y]t,j = sgn(e^_ . - r)efe/'i^'^-, (3.11) 
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where sgn(ji:) = 1 for x ^ and sgn(x) = -1 for jc < 0. By i3.9i and noting that QZ is 
independent of {, we have 

[Sbo-y]ej f,y _f,y 

— — = ecBnZ = ScB„Z ", 

for all 1 ^ (,n ^ r, which is equivalent to saying that P has compatible symmetry. ■ 

Now, for a pair of biorthogonal d-band low-pass filters (ao, ao) with multiplicity 
r satisfying (13.61 ), we have an algorithm (see Algorithm 2) to construct high-pass fil- 
ters ai, . . . , ad-i and ai, . . . , ad-i such that the polyphase matrices P(z) and P(z) de- 
fined as in ( 13.5b satisfy P(z)P*(z) - I^r- Here, Pao [ao;o, • . ■ , ao;d-i] and Pj^ :- 
[ao;o, . ■ . , ao;d-i] the polyphase vectors of ao, ao obtained by (13.2b , respectively. 

Algorithm 2 Construction of Biorthogonal Multiwavelets with Symmetry 

(a) Input: (ao, ao), a pair of biorthogonal d-band filters with multiplicity r and with 
the same symmetry as in (13.6b . 

(b) Initialization: Construct a pair of biorthogonal matrices (U, U) in F by Lemma[3] 
such that both P := PaoU and P = Pao^ (U = (U*)"') are matrices of Laurent 
polynomials with coefficient in F having compatible symmetry: SP - SP - 
{eiz'^^ . . . ,erz'^'Y SO for some ki,...,kr e Z and some 1 x dr row vector 6 of 
Laurent polynomials with symmetry. 

(c) Extension: Derive P^, P^ with all the properties as in Theorem [T] from P, P by 
Algorithmic] _ 

(d) High-pass Filters: Let P := PgU* -: {am;y)o^m,y^d-u P := P^U* =: 
(am;y)o^m,y^d-i as in ( 13.5b . For OT = 1, . . . , d - 1, define high-pass filters 

, d-1 , d-1 

ajz) := - y a„r,y(z'^)z\ a™(z) := - ya,„;y(z'^)z' . (3.12) 

(e) Output: a pair of biorthogonal filter banks ({ao, ai, . . . , aj-i}, {ao, ai, . . . , ad-i)) 
with symmetry and with the perfect reconstruction property, i.e. P, P in (13.5b are 
biorthogonal and all filters am, a,„, m = 1, . . . , d - 1, have symmetry: 



a,„(z) = diag(4'z^^'.", . . . , e';'z^^")am(l/z)diag(fiiz-^' , . . . , SrZ-'''), 
a„,(z) = diag(e'/'zd^'.", . . . , e';zd^")a„(l/z)diag(£iz-^' , . . . , SrZ-''), 



(3.13) 



where c^' := (k^ - k() + C( e R and all e {-1, 1), ;t™ e Z, for f = 1, . . . , r and 
m = 1, . . . , d - 1, are determined by the symmetry pattern of P^ as follows: 

[eiz*' , . . . , srz"'-, slz"' , . . . , slz"' z*?" , . . . , fSe := SP.. (3. 14) 



Let ((p, 4>) be a pair of biorthogonal d-refinable function vectors in L2(K.) associ- 
ated with a pair of biorthogonal d-band filters (ao, ao) and with (f) — . . . , 0^]^, 
(p - [01, ... , (prY ■ Define multiwavelet function vectors i/''" - [i/'™, . . . , il/'yY , i/'™ - 
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[t//", . . ., ij/'"]^ associated with the high-pass filters a,„, a,„, m = 1, . . . , d - 1, by 

r (d^) := a„,(e-'^)m, r (d^) a„(e-'^)0(a ^ e (3.15) 



It is well known that {i/^', 
wavelet basis in LaCI^)- 

Since the high-pass filters ai, . . . , ad-i, ai, . 
ify that each i//"' = [i/r™, . . . , i/r';]^, f" = [0^™, . 
following symmetry: 

^'«(c'« - •) = ^?',>^^ ^2 (^2 - ■) = e^>2 ' 



generates a biorthonormal multi- 

. , acj-i satisfy ( 13.131 ). it is easy to ver- 
. , i/r"']^ defined in ( 13.151 ) also has the 



In the following, let us present several examples to demonstrate our results and 
illustrate our algorithms. 

Example 2. Let d - r = 2 and 00,0^0 be a pair of dual d-filters with symbols ao(z), ao(z) 
(cf. Ill]) given by 



ao(z) = 



16 



^'^'^ ^ 38i 



8 6z-^+6 

-28z"' -H216-28Z 112z-1h-112 
21z-i - 18-H330Z- 18z-H-21z3 -sgr' -h 60 h- 60z - 362^ 



Both ao and ao have the same symmetry pattern and satisfy ( |3.6l l. Let d - d\d2 with 
t/i = I and ii2 = 2. Then, Pag 



p -1 

- 16 

ao 



8 6 
3-z 8 
216 



[ao;o,ao;i] and [a();o,ao;i] are as follows: 

6z-i 



+ 3 
112 



-28(z-i + 1) 



112z- 



-18(1 +z) 12(5 -3z) 3(7z-i -H 110 + 7Z) 12(5 -3z-') 



Let U and U be defined by 



U : = 



1 














1 





1 








1 
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-z _ 



2 

10 1 

2 

z -z 



Then we have UU* = h- Let P := Pa„U and P := Pa„U. Then we have SP = SP 
[l,z]^[l, l,z"', -1] and P, P are given as follows: 



P = 



P = 



4 6 

1(1 -Hz) 4 2(1 -z) 

216 112 -28(1 -I- z"') 

192 [ -18(1 -Hz) 12(l+z) 3(7z-i + 110-)-7z) 48(l-z) 
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Now applying Algorithm 2, we obtain two extension matrices and P^, as follows: 

96 144 

1 24(1 +z) 96 48(1 -z) 
192 _ii2 -3(2-1- 70 + z) -12(1 -6(^1 -z) 

-6(z-z-i) -24(1 -z-i) 12(z+14 + z-') 

216 112 -28(1 +z"') 
1 -18(1 +z) 12(1 +z) 3(7z"' + 110 + 7z) 48(1 -z) 

192 -144 96 -24(1 +z"') 

-96(1 -z"') 192 

Note that SP^ - SPe = [l,z, 1, -1]''^[1, l,z"', -1]. Now from the polyphase matrices 
P :- PfU* -: (a,„;y)o^m,y^i and P :- P^U* -: (am-y)o^m,y^i, we derive two high-pass 
filters ai, ai as follows: 



P. = 



ai(z) 



384 



a,(z) = - 



See Figure 3.1 for the graphs of 

0^= [<Al,«A2]^ 



-8(3z + 28 + 3z">) 3(z^ - 3z + 70 + 70z-' - 3z-^ + z'^) 
-48(z - z"') 6(z^ - 3z + 28 - 28z-' + 3z-^ - z"^) 
-(z + 6 + z-') 4(1 +z-') 
-4(z-z-') 8(1 -z-') 

■■ [(puc/>2f, iff = [lAi.iAa]^, (/> = [(/>\,(f>2V, and 





i \ 



Figure 3.1: The graphs of ip = [if>i,<p2] > lA = [^Ai . 'I'lV (top, left to right), and = [i^i, 1^2] , lA = [i/'i,i)'2Y 
(bottom, left to right) in Example|2] 
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Example 3. Let d = 3, r = 2, and ao, «() be a pair of dual d-filters with symbols 
ao(z), ao(z) (cf. [12]) given by 



a„(z) = ^ 



flii(z) aniz) 
fl2i(z) aiiiz) 



ao(z) 



1 



34884 



where 



flii(z) 
aniz) 
aiiiz) 
aiiiz) 



-212-2 + 30^-1 +81 + 14z- 



an(z) 
aiiiz) 

-5z' 



aniz) 
aiiiz) 



60 z"' +84 



4z 



-2 



•4z-i +84z- 



60 



-5z"' + 14 + 81z + 30z2 



and 



aiiiz) 
aniz) 
aiiiz) 



1292z"2 + 2844 z"' + 17496 + 2590 z ■ 



-4773z"2 + 9682z"' 
-969z"2 + 386z"' 



■21z^ 



1284z2 + 1866z^ 



8715 - 2961 Z + 386Z -969 z' 
■2961 +8715 z + 9682 z^- 4773 z^ 



a22(z) = 1866z"2- 1284 z"' +2590+ 17496 z + 2844 z2 + 1292 zl 
The low-pass filters oq and oo do not satisfy (I3.6l l. However, we can employ a very 
simple orthogonal transform £ := |^ | _', j to ao, ao so that the symmetry in (13.61 1 

holds. That is, for bo(z) := £'ao(z)£'"' and bo(z) '■- £'"'ao(z)£', it is easy to verify that 
bo and bo satisfy ( 13. 6t with ci - C2 = 1/2 and ei = 1,S2 = -1. Let d - d\d2 with 
= 1 and d2 - 3. Construct Pb„ :- [bo;o, bo;i, bo;2] and P^^ := [bo;o, bo;i, bo;2] from bo 
and bo. Let U be given by 



U = 



,-1 





,-1 

^ 








' 


















1 





-1 














1 





-1 




















1 




















1 



and define U (U*) '. Let P PboU and P P^^U. Then we have 5P = 5P = 
[z"', -z"']^[l, -1, -1, 1, 1,-1] and P, P are given by 



P = c 



f2l(l - \) 

T2i(l - i) 



f22(l + \) 
^2(1-7) 
T22(l + i) 



f23(l + \) 
^23(1 + i) 



Ha 

^24(1 - i) 



fl5(l + ^) 
f25(l - 7) 
^5(1 + 7) 
^25(1 - i) 



fl6(l - 7) 
f26(l + 7) 
Ml - 7) 
^26(1 + 7) 



where c = 4^ > c = ™d f^^'s, f^^'s are constants defined as follows: 



486' 

fii = 162; 
t2i = -126; 
Tn = 5814; 
T21 = -5551; 



tl2 
122 

T22 



34884 
34; 
-14; 
-1615; 
5808; 



tl3 
hi 

hi 



-196; 
176; 
-7160; 
7740; 



tl4 
?24 
h4 
T24 



■ 0; 

: -36; 

■ 0; 

: -1358; 



tl5 
hs 

hi 



81; 
-99; 
5814; 
-6712; 



/16 = 29; 
h6 = -31; 
hf, = 2584; 
7,6 = -4254. 
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Applying Algorithm 2, we obtain and as follows: 



?i5(l + i) he{l-k) 



f2l(l-^) f22(l + ^) f23(l+i) r24(l-i) r25(l 



-) f26(l+i) 



f3l(l + ^) f32(l-^) f33(l-j) ?34(l + t) f35(l + j) ?36(1 - j) 



?41 











?44 



^45 











t52 



t53 











^56 



?6l(l-;) + + ?64(l-i) ^65(1-1) ?66(l + i) 

where all ^j^'s are constants given by: 





= 24; 




472 _ 




148 


hi 


^32 - 




^33 






= -36; 


t35 = 


-24; 


t36 


111 




109998 




94041 




109989 


Ul 




t4A = 




Us 




533 ' 


533 ' 


533 ' 


t52 


= 406c(); 


f53 = 


323co; 


he 


= 1 142co; 


tei 


= 24210ci; 


f62 - 


14318ci; 


t(,3 


= -11807ci; 




= -14616ci; 


t66 = 


-1934ci; 


Cl 


= 200/26163 



Co 



1609537 
13122 ' 
-26721ci; 



And 



fll(l + 7) fl2(l- 



Ha 



?15(l + 7) fl6(l-7) 



T2l(l-i) 1^2(1 + ^) ^23(1 + 7) ^(1-7) ^25(1- 7) %(l + 7) 



f3i(i + \) 'h2(\-\) 'hi(\-\) r34(i + j) V35(i + ^) h6{\-\) 



tA5 



hi t53 

where all tj^s are constants given by: 



_ _ _ t56 

■) ?64(l-i) ?65(l-i) f66(l + i) 



f31 


= 3483Fo; 


f32 


= 37427?o; 


^33 = 


4342co; 


?34 = 


-12222co; 


?35 


= -3483?o; 


h6 


= -7267; 


Co = 


8721 

4264' 








= 5814; 


tu 


= 11628; 


t45 = 


-11628; 






hi 


= 3Fi; 


h3 


= 2?i; 


tse = 


lOci; 


Cl = 


12680011 
243 ' 


tei = 


= 18203C2; 


t62 


= 101595C2; 


t63 


= 1638C2; 


t(A 


= -33950C2 


tb5 - 


-- -IO822F2; 


fee 


= -36582?2; 


'ci 


26163 
213200 







Note that Pg and Pg satisfy 

SPe = SPe = \.Z'\ -Z-\z-\ 1, -1, -Z-'Yih -1, -1, 1, 1, -!]• 
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From the polyphase matrices P := PgU* and P 
bi, b2 and bi, as follows: 



bi(z) = 



,b2(z) = 



PgU*, we derive high-pass filters 



b\,{z) b\^{z) 
bUz) bl,(z) 



where 



, 199 



125 , 4 2 

z z + 

6561 81 



199 



125 



6561^ 81^ '*'6561 



-2. 



6561 



b'niz) = 



679 
3198 

bUz) = 



z + 



2132 



z - 



125 

656T 
679 

3198 
387 

2132 



z - 



56 

656T 
679 



z + 



361 



-z + 



6561 6561 



56 125 _2 

z + ^TTT-z ; 



6561 



z - 



z; 



1599 



z ; 



= c3(323z^-323z); 
^z^b) = C3(406z^ + 2284z2 + 406z); 

= C4(-36017+ 12403 -29232 2^ + 36017 z + 29232 z"' - 12403 z"^); 

bj^iz) = C4(41039 - 12403 z^ - 3868 z^ + 41039 z - 3868 z"' - 12403 z"^); 
27 50 



C3 = 



3219074^ 



C4 



6357609' 



And 
where 

bliiz) 
b\xiz) 



bi(z) = 



b\,(z) b\2(z) 
b\,(z) b\,(z) 



,b2(z) = 



bj,(z) bj^iz) 
bl,(z) bl,(z) 



859 7825 

" 17056 17056 ' 

49649 25205 

" 17056 17056 ' 



3483 
8528 



859 

"17056' 
49649 



3483 
8528' 



7825 
17056 ' 



-2. 



559 2 

656^ 17056 



559 

^^ 656^" 



25205 
17056' 



-2. 



i(z'+z-2z2); bl^iz) ^ W - z); 
6 3 



^ti(z) = 2c3(z'-z); 



^^2(z)=?-3(3z^ + 10z' + 3z);^3 = 



39257 



^21 W = - 



9939 



+ 



59523 



^22(^) 



170560 852800 
81327 40587 



z^ - 



170560 170560 



z'- 



26244' 
16233 

426400 

4221 , 
z + 



z" + 



9939 



■z + 



16233 



59523 



170560 426400 



z - 



852800 



-2. 



32800 



81327 
170560 



z - 



4221 
32800 ' 



-1 



40587 



170560 



Then the high-pass filters bi,b2 and bi,b2 satisfy ( 13.131 ) with cj 



1/2, si 



3/2, e[ = -1,62 ~ respectively. 



1 , fij = 1 ™d c J = c: 

Let 31,32 and 3i,32 be high-pass filters constructed from bi,b2 and bi,b2 by 
3i(z) := £-'bi(z)£,32 := £-'b2£ and 3i(z) £bi(z)£"\'32 Eb2E-\ 
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See Figure 3.2 for the graphs of the 3-refinable function vectors 0,0 associated 
with the low-pass filters ao, ao, respectively, and the biorthogonal multiwavelet func- 
tion vectors 0\0^ and i/f\0^ associated with the high-pass filters a 1,82 and a 1,82, 
respectively. Also, see Figure 3.3 for the graphs of the 3-refinable function vectors 
77, 77 associated with the low-pass filters bo, bo, respectively, and the biorthogonal mul- 
tiwavelet function vectors and associated with the high-pass filters bi, b2 
and bi, b2, respectively. 




Figure 3.2: The graphs of = [rn,ri2f, = and f = [({,(1^ (top, left to right), and^ = 

= aHlY, and?^ = [fj,:^]^ (bottom, left to right). 




Figure 3.3: The graphs of = [0i,02]^, iA' = , and ijj^ = (top, left to right), and 

J= 0u^2f,^^ = [ifrj.l^'j]^, and(/^^ = l^j.^lV (bottom, left to right). 
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4. Conclusions and Remarks 



In this paper, we study the matrix extension problem with symmetry for the bio- 
thogonal case. We obtain a result on representing a pair ofrxs biorthogonal matrices 
(P, P) having the same compatibly symmetry and provide a step-by-step algorithm for 
deriving a pair of s x s biorthogonal matrices from a given pair of biorthogonal matri- 
ces (P, P). Our results show that for the one row case (r = 1), the support lengths of 
the extension matrices can be controlled by the given pair of columns. We apply our 
results in this paper to the derivation of synnmetric biortahogonal multiwavelets from a 
pair of dual d-refinable functions. 
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